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The  rate  of  industrial  accident  occurrence  has  remained  a  perennial 
challenge.  This  paper  points  to  the  need  for  the  application  of  Markov 
chain  model  in  unwrapping  the  deeper  meanings  buried  in  the  safety 
incident  data  that  mere  descriptive  statistics  can  hardly  furnish.  In  line 
with  the  study  design,  a  16-year  data  was  obtained  from  Nigeria  Gas 
Company,  a  subsidiary  of  Nigeria  National  Petroleum  Corporation 
(NNPC).  The  historical  incidence  records  were  characterized  and 
proved  to  possess  a  note  of  stochastic  regularity  that  fits  into  a  Markov 
Chain  Model.  A  twenty  state  transition  was  used  for  the  study,  namely: 
fatality,  third  party  fatality,  permanent  disability  for  example.  Result 
emanating  from  the  study  reveals  that  subjects  make  about  thirteen 
habituations  among  various  states  in  the  organization  before  being 
absorbed  in  any  ten  absorbing  states  with  a  standard  deviation  of  12. 
Remarkably,  70.5%  of  the  field  workers  in  the  organization  had 
noteworthy  severe  medical  treatment  case.  In  conclusion,  the  Markov 
Chain  Model  was  able  to  identify  states  such  as  unsafe  acts  and  unsafe 
conditions  transitions  to  have  influenced  incident  levels  the  most  in  the 
organization.  This  study  has  also  shown  that  Markov  Chain  model  can 
be  successfully  applied  to  industrial  accident  data,  unveiled  significant 
visits,  habituations  which  the  organization  can  explore  in  optimizing 
their  injury  prevention  programme  and  ensures  field  staff  safety. 


1.  Introduction 

In  spite  of  the  multi-strategic  approach  to  nipping  industrial  accident  occurrences  to  the  base  in 
most  factories,  the  phenomenon  appears  progressively  unabating  of  which  the  oil  and  gas 
industry  situation  is  no  exception.  This  is  attributed  to  manager’s  lack  of  commitment  and 
dependence  on  mere  descriptive  statistics  as  a  major  tool  for  analyzing  and  monitoring  accident 
statistical  data  to  establish  trend  and  patterns.  Moreover,  [1]  corroborated  this  fact,  according  to 
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them  only  with  descriptive  statistics  it  is  not  possible  to  carry  out  an  in  depth  review  of  the 
causes  of  accidents.  Therefore,  statistical  methodologies  have  been  improving  which  have 
enabled  better  safety  design  and  policy  improvement.  Hence,  the  need  for  new  ways  and  better 
robust,  rigorous  engineering  models  such  as  the  Markov  chain  model  which  is  considered 
effective  for  analyzing  accident  statistical  historic  records  in  order  to  unveil  deeper  meanings 
embedded  in  such  accident  data  aimed  at  reducing  accident  risk  occurrence  that  descriptive 
statistics  can  hardly  provide.  The  Markov  chain  model  has  the  capability  of  not  just  establishing 
trend  and  patterns  but  unwrapping  insightful  transitions,  habituation  workers  in  an  organization 
undergo,  movements,  visit  of  certain  state  events  before  they  are  absorbed  by  a  major  life 
threatening  state. 

However,  accident  is  yet  under  reported,  arbitrary,  and  in  some  cases  in  Nigeria  highly 
customized.  According  to  [2]  fatality,  injury  and  accident  rates  are  increasing  and  management 
commitments  to  accident  remain  poor  and  this  arguably  contributes  largely  to  accident  causation 
accounting  for  91.3%  contributory  factor.  It  has  been  observed,  over  the  years,  that  manager’s  of 
oil  and  gas  industries  are  more  reactive  to  health  and  safety  issue  as  they  occur  than  before  they 
are  likely  to  happen.  The  foregoing  became  obvious  when  investigations  into  several  episodic 
accidents  such  as  the  Piper  Alpha  disaster  (1988),  Bhopal  Gas  plant  disaster  (1984),  Chernobyl 
Nuclear  power  plant  disaster  (1986)  and  Deep  Water  Horizon  Oil  Spill  disaster  (2010) 
highlighted  the  role  of  management  negligence  and  the  relevance  of  organizational  level,  lack  of 
proper  knowledge  and  insincerity  shown  as  pivotal  to  these  catastrophic  outcomes  [3,  4,  5], 
Therefore  developing  better  proactive  approach  and  models  to  spotlight  dangerous  occurrences 
through  a  veritable  decision  supportive  model  becomes  necessary.  As  posited  by  [6],  if 
inadequate  intervention  patterns  are  adopted  accidents  proliferate. 

The  Markov  chain  models  have  been  successfully  applied  in  other  fields  of  research  work  in 
engineering.  As  asserted  by  [7]  Markov  method  can  be  powerful  tool  in  reliability, 
maintainability  and  safety  (RMS)  engineering.  Markov  chains  are  commonly  applies  to  the  study 
of  dependability  of  complex  system.  While  [8,9],  avers  to  the  fact  that  Markov  models  are 
powerful  statistical  tool  and  do  have  a  joint  history  and  form  a  fruitful  partnership  with 
maintenance  modeling,  also  been  successfully  applied  in  component  diagnostics,  prognostics  and 
maintenance  optimization  across  a  gamut  of  industries. 

Furthermore,  the  applications  of  this  model  in  the  area  of  human  life  safety  in  the  oil  and  gas 
sector  are  limited  and  therefore  need  further  attention.  Markov  chain  model  was  applied  by  [10] 
in  Robot  safety  identifying  potential  risk  for  industrial  robot  and  the  definition  of  hazard  rate  at 
different  state  for  robot  system.  Again  [11]  reports  a  fruitful  application  of  Markov  chain  in 
predicting  risk  severity  and  exposure  level  of  workers  in  Warri  Refining  and  Petrochemical 
Company  (WRPC)  involving  four  states,  with  two  absorbing  states  and  two  non-absorbing  states. 
A  survey  approach  involving  that  of  questionnaire  administration  instrument  to  fifty  workers  was 
adopted.  Again,  [12]  did  a  longitudinal  study  appraising  the  patterning  of  episodic  incidence  of 
industrial  accidents  in  oil  and  gas  firm  in  the  Niger-delta  area  of  Nigeria  using  10-years  historical 
data.  Essentially,  this  study  widens  the  horizon  of  the  application  of  the  Markov  chain  model 
covering  a  20-State  system,  with  ten  absorbing  and  non-absorbing  states  respectively. 

In  general,  the  study  is  aimed  at  unwrapping  the  deeper  meanings,  important  prediction  or 
implications  buried  in  the  industrial  accident  records  by  applying  Markov-chain  model  as  a 
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veritable  tool  that  will  aid  manager’s  decision  in  policy  making  and  improving  health,  safety  and 
environmental  performances. 

2.  Methodology 

In  an  attempt  to  decipher  the  deeper  meanings  buried  in  the  industrial  accident  data  of  Nigerian 
Gas  Company  Limited  (NGC),  a  subsidiary  of  Nigerian  National  Petroleum  Company  (NNPC), 
and  a  Markov  chain  molded  was  fitted  into  a  16-year  (2000  -  2015)  historic  safety  data  as 
depicted  in  Table  1.  The  industrial  accident  data  were  examined  for  embedded  Markov 
properties  namely  stochastic  regularity,  absorbing  behaviour  and  the  long-run  distribution 
amongst  the  various  states.  The  basic  assumptions,  applicable  theorems  and  formulae  leading  to 
the  computations  and  decision  making  process  are  precisely  stated.  Furthermore,  the  statistical 
computations  were  done  with  the  aid  of  (MATLAB  .R.  2016  a)  software  from  which  deduction 
and  inference  were  derived  from  the  results  obtained  and  which  guided  subsequent  discussion 
and  conclusion. 


Table  1.  NGC  HSE  Historical  Record  from  NNPC  (2000  -  2015) 


S/N 

States 

Total 

1. 

Fatality  (FT) 

- 

Absorbing  States 

2. 

Third  Party  Fatality  (TFF) 

1 

3. 

Permanent  Disability  (PD) 

1 

4. 

Fire  Incidence  (severe)  (FIs) 

23 

5. 

Injury  (severe)  (Is) 

13 

6. 

Medical  Treatment  Case  (severe)  MTCs 

165 

7. 

Road  Traffic  Accident  (severe)  RTAs 

- 

8. 

Oil  Spill  (callouts)  OS 

- 

9. 

Restricted  Work  Case  (severe)  RWCs 

6 

10. 

Human  Error  (HE) 

16 

11. 

Lost  Time  Days  LTD 

30 

i 

Ion-Absorbing  States 

12. 

Injury  (minor)Im 

105 

13. 

Fire  Incidence  (minor)  FIm 

101 

14. 

Medical  Treatment  Case  (minor)  MTCm 

1202 

15. 

Restricted  work  case  (minor)  RWCm 

11 

16. 

Road  Traffic  Accident  (minor)  RTAm 

9 

17. 

First  Aid  Case  (FAC) 

32 

18. 

Unsafe  Condition  (UC) 

365 

19. 

Unsafe  Acts  (UA) 

368 

20. 

Near  Misses  (NM) 

104 

Total 

2552 

^ ~ 

2.1.  Theoretical  Formulation 
(a)  Lemmatization 

In  mathematical  theory  of  probability,  a  Markov  chain  is  an  absorbing  chain  provided: 
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Lemma  1:  there  are  one  or  more  absorbing  states.  Furthermore,  a  state  in  a  Markov  chain  is 
said  to  be  absorbing  if  the  probability  of  an  object  leaving  the  state  once  entered  is  zero,  or 
once  entered  cannot  be  exited  and  the  probability  that  it  stays  in  that  state  is  one. 

Lemma  2:  from  each  of  the  non-absorbing  states,  it  is  possible  to  reach  some  absorbing  state 
in  the  long  run  by  a  number  of  steps.  In  an  absorbing  Markov  chain,  a  state  that  is  not 
absorbing  i.e  non-absorbing  is  also  called  transient. 

(b)Theorem 

In  standard  form,  the  transition  matrix  T  is  expressed  as 


T  = 


Abs 

Non-abs 


Abs  Non-abs 


(I 

0\ 

u 

QJ 

(1) 


Then  in  standard  form: 

a.  The  matrix  (I-Q)  is  invertible. 

Where  1=  identity  matrix  and  Q  is  non-absorbing  matrix.  In  other  words  Q  is  a 

Matrix  having  an  initial  distribution  (or  describe  the  probability  of  transitioning  from  some 

transient  state  (non-absorbing  state)  to  another. 

0  =  zero  matrix  or  null  matrix,  while  R  matrix  describes  the  probability  of  transitioning  from 
some  non-absorbing  state  (transients  state)  to  some  absorbing  state. 

For  the  purpose  of  this  study,  we  seek  to  state  the  long-run  distribution  of  T  i.e  Tn  =  T  from 
which  the  fundamental  matrix  (N)  is  derived  without  proof. 


rj^n  _ 


0 


[l  +  Q+Q1  +Q"  +...+Q"  ')R  W2 


(2) 


The  long  run  distribution  of  T 

r=f=(.  '  ;°) 

((/-eT'tfio)  (3) 

Equation  (iii)  is  key  to  the  computational  analysis  that  follows,  and  from  which  the  fundamental 
matrix  N  =  (I  -  Q)'1  is  derived 

(c)  Fundamental  Matrix  (N), 

N  =  (I  -  Q)'1  (4) 

This  is  a  basic  property  of  an  absorbing  Markov  chain  and  it  represents  the  expected  or  average 
number  of  visits  or  number  of  times  objects  starts  in  the  ith  non- absorbing  state  to  a 
state(transient)  j  ’  before  being  absorbed. 

Primarily,  with  fundamental  matrix  (N)  in  hand,  other  properties  or  derivatives  are  obtained, 
such  as  variance  on  number  of  visits  N2=  N  (2Ndg  -  I)  -  Nsq,  expected  number  of  step  or 
cumulative  movement  of  object  r  =  N  variance  on  number  of  steps 
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T2=  (2N  -  I)  t  -  rsq,  the  long  run  distribution  of  object  or  staff  among  the  various  absorbing  states 
B  =  NR  and  the  transient  probabilities  H=  (N  - 1)  Ndg_1 


(d)  Variance  on  number  of  visits  (N2), 


N2  =  N(2Ndg-I)-Nsq  (5) 

This  gives  the  associated  variance  on  the  expected  number  of  visits  or  simply  the  variance  on  the 
number  of  visits  to  a  non-absorbing  state  (transient)  j  starting  at  a  transient  state  i  (before  being 
absorbed)  is  the  (ij)-  entry  of  the  matrix. 

where  in  Equation  (5),  N  is  the  fundamental  matrix  (N),  Ndg  is  diagonal  of  N  with  zeroes  and  Nsq 
is  the  square  of  each  element  of  the  N  matrix. 

(e)  Expected  number  of  steps  or  movements  (t) 


t  =  N  ^ 


where  = 


(1\ 

1 

1 

1 

1 

1 

1 

1 

1 


is  a  column  vector  whose  entries  are  all  unity 


(6) 


The  expected  number  of  steps  or  total  number  of  movements  of  staff  or  objects  within  the  non¬ 
absorbing  state  before  being  absorbed. 

(f)  Variance  on  number  of  step  or  movements. 

T2=(2N-I)T-Tsq  (7) 

The  variance  on  the  number  of  steps  or  movement  before  being  absorbed  when  starting  in  a  non¬ 
absorbing  state  i.  or  termed  the  associated  variance,  where  tsq  is  the  hadamard  product  of  t  with 
itself  (i.e.  each  entry  of  t  is  squared). 

(g)  Absorbing  Probabilities  ( B ) 

B  =  (  I-Q)1R  =  NR  (8) 

This  specifies  the  long-run  distribution  of  objects  among  the  various  absorbing  state,  provided  all 
object/staff  start  in  the  non-absorbing  state. 

(h)  Transient  Probabilities  (H) 

H=  (N  -  I)  Ndg_1  (9) 

This  gives  the  probability  of  likelihood  or  chances  of  visiting  non-absorbing  state  j  starting  at  a 
non-absorbing  state  i. 

2.2.  Statistical  Computations 

The  statistical  computation  when  Markov  chain  is  applied  involves  two  (2)  regimes  viz: 

I.  the  absorbing  state  regime;  and 

II.  the  non-absorbing  state  regime. 
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2.2.1.  Absorbing  State  Computations 

The  diagraph  of  the  20-states  structure  in  Fig  1  gives  an  insight  on  the  way  the  system  operates. 
In  this  regime  all  the  transition  probabilities  associated  with  the  absorbing  states  are  heuristically 
determined.  Basically,  heuristic  algorithm  method  is  adopted  which  involves  reasoning  in  the 
determination  of  the  absorbing  state  probabilities.  For  example  Pn  =  Pr(FT)  =  1 
Mathematically  this  can  be  expressed  as  Pn  =  Pr(FT)  =  1,  implying  that  a  staff  who  is  wasted  by 
fatality  (death)  remains  dead.  This  is  a  case  of  persistence.  As  stated  above,  it  follows  similar 
heuristic  determination  applies  for  others  and  can  be  expressed  mathematically  as  P22  =  P33  =  P44 

. =  P10, 10  =  1-  Furthermore,  P12  which  represents  a  subject  transiting  from  fatality  (FT)  to  third 

party  fatality  (TPF);  it  is  a  case  which  is  implausible;  in  other  words,  this  simply  means 
transiting  from  a  state  considered  absorbing  to  another  state  which  is  an  absorbing  or  non¬ 
absorbing  is  implausible.  Then,  similar  heuristic  argument  are  adopted,  P12  =  P21  =  P13  =  P31  = 
P 14  =  P41  =  . . .  =  Pi, 10  =  Pkm  for  reasons  of  Implausibility. 


Connectors  from 
absorbing  to 
non-absorbing  states 


Absorbing 

Non-absorbing 


Fig.l.  Diagraph  of  the  20-State  structure  each  state  is  linked  to  all  the  20-states  by  an 
outflow  and  return  lines 
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2.2.2  Non-absorbing  State  Computation 

In  the  determination  of  probabilities  under  the  non-absorbing  state  regime,  Bayesian  method  is 
adopted.  For  this  study,  ten  (10)  non-absorbing  states  are  considered  namely,  Lost  Time  Days 
(LTD),  Injury  (minor)  Im,  Fire  incidence  (minor)  FIm,  Medical  treatment  case  (minor)  MTCm, 
Restricted  work  case  (minor)  RWCm,  Road  traffic  accident  (minor)  RTAm,  FAC,  UC,  UA,  and 
NM  as  shown  in  Table  1.  In  the  use  of  Bayesian  methodology  to  compute  transition  probabilities 
under  non-absorbing  regime  certain  assumptions  were  considered  to  allow  for  mathematic 
tractability.  The  assumptions  follow. 

Assumption  1:  A  staff  that  is  considered  by  management  to  have  a  lost  time  day  cannot  at  the 
same  time  be  available  to  commit  near  miss  (NM=  0),  Pn,  20- 

Assumption  2:  A  staff  who  had  just  suffered  a  minor  injury  (Im)  as  an  entry,  cannot  be  said  to 
have  minor  fire  incidence  at  the  same  time  (FIm  =0),  Pi2,i3=0 

Assumption  3:  A  well  trained  staff  that  had  a  minor  fire  incidence  (FIm)  in  the  industry  as 
recorded  cannot  be  involved  in  a  minor  road  traffic  accident  at  the  same  time  (a  case  of 
impossibility)  RTAm=0,  P 1 3,1 6=0- 

Assumption  4:  A  staff  whose  state  is  considered  as  a  minor  medical  treatment  case  MTCm 
cannot  be  said  to  be  in  a  state  of  permanent  disability  PD  =0,  Pm,  3=0. 

Assumption  5:  A  staff  who  at  a  particular  time  is  in  a  state  considered  as  a  minor  restricted  work 
case  (RWCm)  will  be  unavailable  to  be  exposed  to  unsafe  condition  (UC  =  0),  P15 18=0. 
Assumption  6:  A  healthy  contract  staff  who  had  a  minor  road  traffic  accident  (RTAM)  cannot 
transit  instantly  to  third  party  fatality  (death)  TPF=0,  Pi6,2=  0. 

Assumption  7:  A  staff  whose  injury  is  considered  to  be  a  first  aid  case  (FAC)  cannot  be  said  to 
have  committed  a  major  work  related  human  error  (HE  =  0),  P1740 

Assumption  8:  It  is  impossible  for  a  staff  to  be  at  a  state  considered  an  unsafe  condition  (UC)  at 
the  same  time  be  said  to  be  in  severe  restricted  work  case  (RWCs  =  0),  Pis, 9=  0 
Assumption  9:  A  well  trained  staff  who  is  involved  in  an  unsafe  act  (UA)  cannot  instantly  transit 
to  a  state  of  severe  medical  treatment  case  MTCs  =  0,  Pi9,6  =  0 

Assumption  10:  A  well  trained  staff  who  commits  near  miss  cannot  instantly  transit  to  fatality 
FT=  0,  P2o,i  =  0 

Basically,  with  the  ten  (10)  non-absorbing  state  considered  previously,  a  ten  different 
combinations  of  the  20  state  structure  were  arranged  in  Table  2  with  the  columns  indicating 
different  combination.  Furthermore,  in  applying  the  Bayesian  approach  it  uses  conditional 
probability  of  events  whose  occurrences  are  presumptuously  equated  to  zero  in  every  column 
except  those  whose  state  total  is  zero  as  shown  in  Table  1  represented  as  serial  number  1,  7,  8 
which  are  fatality,  road  traffic  accident  severe  and  oil  spill  respectively. 

Column  I(Vseti)  for  instance,  the  probability  of  subject  who  commits  Unsafe  Act  (UA) 
transiting  to  First  Aid  Case  (FAC),  given  that  the  condition  that  subject  has  not  had  near  miss 
earlier,  denoted  mathematically  as  P{(UA  — ►  (FAC)lNM  =  0},  similarly,  in  the  same  column  the 
probability  of  subject  who  is  exposed  to  unsafe  condition  transitioning  to  the  commission  of 
human  error,  on  the  account  that  the  subject  has  not  had  near  miss  earlier,  is  stated  thus 
P{(FAC— >HElNM=  0)}. Similar  representations  apply  for  all  cases  under  column  l(V-set  1). 
Computations  for  all  the  sets  are  collated  and  depicted  in  Table  3  also,  see  Table  4  for  mode  of 
transitions 
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Table  2.  Ten  Different  Combination  of  the  20  state  structure  (NGC)  are  arrange  as  follows 


S/NO 

V_SET  1 

V_SET  2 

V_SET  3 

V_SET  4 

V_SET  5 

V_SET  6 

V_SET  7 

V_SET  8 

V_SET  9 

V_SET  10 

1 

FT=0 

11=0 

11=0 

FT=0 

11=0 

FT=0 

11=0 

11=0 

FT=0 

11=0 

2 

TPF=1 

TPF=1 

TPF=1 

TPF=1 

TPF=1 

TPF=0 

TPF=1 

TPF=1 

TPF=1 

TPF=1 

3 

PD=1 

PD=1 

PD=1 

PD=0 

PD=1 

PD=1 

PD=1 

PD=1 

PD=1 

PD=1 

4 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

FIs=23 

5 

Is=13 

Is=13 

Is=  1 3 

Is=13 

Is=  1 3 

Is=13 

Is=  1 3 

Is=  1 3 

Is=13 

Is=  1 3 

6 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=165 

MTCs=0 

MTCs=165 

7 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

RTAs=0 

8 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

os_co=o 

9 

RWCs=6 

RWCs=6 

RWCs=6 

RWCs=6 

RWCs=6 

RWCs=6 

RWCs=6 

RWCs=0 

RWCs=6 

RWCs=6 

10 

HE=16 

HE=16 

HE=16 

HE=16 

HE=16 

HE=16 

HE=0 

HE=16 

HE=16 

HE=16 

11 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

LTD=30 

12 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

Im=105 

13 

FIm=101 

Flm=0 

FIm=101 

FIm=101 

FIm=101 

FIm=101 

FIm=101 

FIm=101 

FIm=101 

FIm=101 

14 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

MTCm=1202 

15 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

RWCm=ll 

16 

RTAm=9 

RTAm=9 

RTAm=0 

RTAm=9 

RTAm=9 

RTAm=9 

RTAm=9 

RTAm=9 

RTAm=9 

RTAm=9 

17 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

FAC=32 

18 

UC=365 

UC=365 

UC=365 

UC=365 

uc=o 

UC=365 

UC=365 

UC=365 

UC=365 

UC=365 

19 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

UA=368 

20 

NM=0 

NM=104 

NM=104 

NM=104 

NM=104 

NM=104 

NM=104 

NM=104 

NM=104 

NM=104 

TOTAL 

2448 

2451 

2543 

2551 

2187 

2551 

2536 

2546 

2387 

2552 

In  column  5(Vset5):  The  probability  of  a  subject  who  had  a  minor  Fire  Incidence  (FIm)  in  row  13 
transiting  to  commission  human  error  HE  in  column  10  (i-j)  entry  on  the  condition  that  he  has 
not  been  exposed  to  an  unsafe  condition,  expressed  mathematically,  P{(FIM— ►  HEjlUC  =  0}. 
Similar  representation  determination  for  each  of  the  columns  is  done  following  similar  pattern 
using  conditional  probability. 

The  Bayesian  approach  essentially  uses  ratio  of  subject  in  any  state  to  the  total  population  of  the 
ten  different  combinations  in  Table  2.  For  instance,  column  Vsetl  serial  number  1  fatality  (FT)  = 
0  which  a  subject  state  1,  to  the  ratio  of  total  population  2448.  This  is  expressed  mathematically 
as  depicted  in  Table3  probability  set  1.  Again,  similar  computations  are  done  for  all  Vsetl  and 
other  Vsets  in  Table  2  and  mathematically  computed  to  obtain  probability  setl  and  other 
probability  sets  in  Table3. 
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Table  3.  Sample  Computation  ofTransition  Probabilities^  using  the  Bayesian  Methodology k 


Probabilities  Setl 

Probabilities  Set2 

Probabilities  Set3 

Probabilities  Set4 

P(FT)=(0/2448)=0.0000 

P(FT)=(0/245 1  )=0.0000 

P(FT)=(0/2543)=0.0000 

P(FT)=(0/255 1  )=0.0000 

P(TPF)=(  l/2448)=0.0004 

P(TPF)=(  1/245 1  )=0.0004 

P(TPF)=(  l/2543)=0.0004 

P(TPF)=(  1/255 1  )=0.0004 

P(PD)=(l/2448)=0.0004 

P(PD)=(  1/245 1)=0.0004 

P(PD)=(l/2543)=0.0004 

P(PD)=(0/255 1)=0.0000 

P(FI_S)=(23/2448)=0.0094 

P(FI_S)=(23/245 1  )=0.0094 

P(FI_S)=(23/2543)=0.0090 

P(FI_S)=(23/255 1)=0.0090 

P(I_S)=(  1 3/2448)=0.0053 

P(I_S)=(  1 3/245 1  )=0.0053 

P(I_S)=(  1 3/2543)=0.005 1 

P(I_S)=(  1 3/255 1  )=0.005 1 

P(MTC_S)=(165/2448)=0.0674 

P(MTC_S)=(  1 65/245 1  )=0.067 3 

P(MTC_S)=(165/2543)=0.0649 

P(MTC_S)=(165/255 1)=0.0647 

P(RTA_S)=(0/2448)=0.0000 

P(RTA_S)=(0/245 1  )=0.0000 

P(RTA_S)=(0/2543)=0.0000 

P(RTA_S)=(0/255 1)=0.0000 

P(OS_CO)=(0/2448)=0.0000 

P(OS_CO)=(0/245 1  )=0.0000 

P(OS_CO)=(0/2543)=0.0000 

P(OS_CO)=(0/255 1)=0.0000 

P(RWC_S)=(6/2448)=0.0025 

P(RWC_S)=(6/245 1)=0.0024 

P(RWC_S)=(6/2543)=0.0024 

P(RWC_S)=(6/255 1)=0.0024 

P(HE)=(16/2448)=0.0065 

P(HE)=(  1 6/245 1  )=0.0065 

P(HE)=(  1 6/2543)=0.0063 

P(HE)=(  1 6/255 1  )=0.0063 

P(LTD)=(30/2448)=0.0123 

P(LTD)=(3 0/245 1  )=0.0 1 22 

P(LTD)=(30/2543)=0.01 18 

P(LTD)=(30/255 1  )=0.0 118 

P(I_M)=(105/2448)=0.0429 

P(I_M)=(  1 05/245 1  )=0.0428 

P(I_M)=(  1 05/2543)=0.04 1 3 

P(I_M)=(  1 05/255 1  )=0.04 1 2 

P(FI_M)=(  1 0 1/2448)=0.04 1 3 

P(FI_M)=(0/245 1  )=0.0000 

P(FI_M)=(  1 0 1/2543)=0.0397 

P(FI_M)=(  101/2551  )=0.0396 

P(MTC_M)=(  1 202/2448)=0.49 1 0 

P(MTC_M)=(  1202/245 1)=0.4904 

P(MTC_M)=(1202/2543)=0.4727 

P(MTC_M)=(  1 202/255 1  )=0.47 1 2 

P(RWC_M)=(1 1/2448)=0.0045 

P(RWC_M)=(1 1/245 1)=0.0045 

P(RWC_M)=(1 1/2543)=0.0043 

P(RWC_M)=(  1 1/255 1  )=0.0043 

P(RTA_M)=(9/2448)=0.0037 

P(RTA_M)=(9/245 1)=0.0037 

P(RTA_M)=(0/2543)=0.0000 

P(RTA_M)=(9/2551)=0.0035 

P(FAC)=(32/2448)=0.0 131 

P(FAC)=(32/245 1  )=0.0 1 3 1 

P(FAC)=(32/2543)=0.0126 

P(FAC)=(32/255 1  )=0.0 1 25 

P(UC)=(365/2448)=0. 1491 

P(UC)=(365/245 1)=0. 1489 

P(UC)=(365/2543)=0. 1435 

P(UC)=(365/255 1)=0. 143 1 

P(UA)=(368/2448)=0. 1 503 

P(UA)=(368/245 1)=0. 1501 

P(UA)=(368/2543)=0. 1447 

P(UA)=(368/255 1)=0. 1443 

P(NM)=(0/2448)=0.0000 

P(NM)=(  1 04/245 1  )=0.0424 

P(NM)=(104/2543)=0.0409 

P(NM)=(  1 04/255 1  )=0.0408 

Table  4.  NGC  mode  of  Transition  Matrix  Probability  Tableau 


FT 

TPF 

PD 

Fis 

Is 

MTCs 

RTAs 

Osco 

RWCs 

HE 

LTD 

Im 

Fim 

MTCm 

RWCm 

RTAm 

FAC 

UC 

UA 

NM 

FT 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TPF 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

PD 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Fis 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Is 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

MTCs 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

RTAs 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Osco 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

RWCs 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

HE 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

LTD 

0 

0.0004 

0.0004 

0.0094 

0.0053 

0.0674 

0 

0 

0.0025 

0.0065 

0.0123 

0.0429 

0.0413 

0.491 

0.0045 

0.0037 

0.0131 

0.1491 

0.1503 

0 

Im 

0 

0.0004 

0.0004 

0.0094 

0.0053 

0.0673 

0 

0 

0.0024 

0.0065 

0.0122 

0.0428 

0 

0.4904 

0.0045 

0.0037 

0.0131 

0.1489 

0.1501 

0.0424 

Fim 

0 

0.0004 

0.0004 

0.009 

0.0051 

0.0649 

0 

0 

0.0024 

0.0063 

0.0118 

0.0413 

0.0397 

0.4727 

0.0043 

0 

0.0126 

0.1435 

0.1447 

0.0409 

MTCm 

0 

0.0004 

0 

0.009 

0.0051 

0.0647 

0 

0 

0.0024 

0.0063 

0.0118 

0.0412 

0.0396 

0.4712 

0.0043 

0.0035 

0.0125 

0.1431 

0.1443 

0.0408 

RWCm 

0 

0.0005 

0.0005 

0.0105 

0.0059 

0.0754 

0 

0 

0.0027 

0.0073 

0.0137 

0.048 

0.0462 

0.5496 

0.005 

0.0041 

0.0146 

0 

0.1683 

0.0476 

RTAm 

0 

0 

0.0004 

0.009 

0.0051 

0.0647 

0 

0 

0.0024 

0.0063 

0.0118 

0.0412 

0.0396 

0.4712 

0.0043 

0.0035 

0.0125 

0.1431 

0.1443 

0.0408 

FAC 

0 

0.0004 

0.0004 

0.0091 

0.0051 

0.0651 

0 

0 

0.0024 

0 

0.0118 

0.0414 

0.0398 

0.474 

0.0043 

0.0035 

0.0126 

0.1439 

0.1451 

0.041 

UC 

0 

0.0004 

0.0004 

0.009 

0.0051 

0.0648 

0 

0 

0 

0.0063 

0.0118 

0.0412 

0.0397 

0.4721 

0.0043 

0.0035 

0.0126 

0.1434 

0.1445 

0.0408 

UA 

0 

0.0004 

0.0004 

0.0096 

0.0054 

0 

0 

0 

0.0025 

0.0067 

0.0126 

0.044 

0.0423 

0.5036 

0.0046 

0.0038 

0.0134 

0.1529 

0.1542 

0.0436 

NM 

0 

0.0004 

0.0004 

0.009 

0.0051 

0.0647 

0 

0 

0.0024 

0.0063 

0.0118 

0.0411 

0.0396 

0.471 

0.0043 

0.0035 

0.0125 

0.143 

0.1442 

0.0408 
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3.  Results  and  Discussion 

Standard  form  of  transition  matrix 


(1 

O') 

u 

QJ 

The  computed  transition  probabilities  which  are  consistent  with  Equation  (1)  and  as  displayed  in 
the  Table  4,  are  depicted  in  the  accompanying  matrix  below.. 

The  Transition  probabilities 


(1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  " 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

"o~ 

0.0004 

0.0004 

0.0094 

0.0053 

0.0674 

T 

~o~ 

0.025 

0.0065 

0.0123 

0.0429 

0.0413 

0.491 

0.0045 

0.0037 

0.0131 

0.1491 

0.1503 

0 

0 

0.0004 

0.0004 

0.0094 

0.0053 

0.0673 

0 

0 

0.024 

0.0065 

0.0122 

0.0428 

0 

0.4904 

0.0045 

0.0037 

0.0131 

0.1489 

0.1501 

0.0424 

0 

0.0004 

0.0004 

0.009 

0.0051 

0.0649 

0 

0 

0.024 

0.0063 

0.0118 

0.0413 

0.0397 

0.49727 

0.0043 

0 

0.0126 

0.1435 

0.1447 

0.0409 

0 

0.0004 

0 

0.009 

0.0051 

0.0647 

0 

0 

0.024 

0.0063 

0.0118 

0.0412 

0.0396 

0.4712 

0.0043 

0.0035 

0.0125 

0.1431 

0.1443 

0.0408 

0 

0.0005 

0.0005 

0.00105 

0.0059 

0.0754 

0 

0 

0.027 

0.0073 

0.0137 

0.048 

0.0462 

0.5496 

0.005 

0.0041 

0.0146 

0 

0.1683 

0.0476 

0 

0 

0.0004 

0.009 

0.0051 

0.0647 

0 

0 

0.024 

0.0063 

0.0118 

0.0412 

0.0396 

0.4712 

0.0043 

0.0035 

0.0125 

0.1431 

0.1443 

0.0408 

0 

0.0004 

0.0004 

0.0091 

0.0051 

0.0651 

0 

0 

0.024 

0 

0.0118 

0.0414 

0.0398 

0.474 

0.0043 

0.0035 

0.0126 

0.1439 

0.1451 

0.041 

0 

0.0004 

0.0004 

0.009 

0.0051 

0.0648 

0 

0 

0 

0.0063 

0.0118 

0.0412 

0.0397 

0.4721 

0.0043 

0.0035 

0.0126 

0.1434 

0.1445 

0.0408 

0 

0.0004 

0.0004 

0.0096 

0.0054 

0 

0 

0 

0.025 

0.0067 

0.0126 

0.044 

0.0423 

0.5036 

0.0046 

0.0038 

0.0134 

0.1529 

0.1542 

0.0436 

lo 

0.0004 

0.0004 

0.009 

0.0051 

0.0647 

0 

0 

0.024 

0.0063 

0.0118 

0.0411 

0.0396 

0.471 

0.0043 

0.0035 

0.0143 

0.143 

0.1442 

0.0408y 

0 


(10) 

Q 


3.1  Computation  of  the  fundamental  matrix 

n=(i-qY 

Computation  of  the  fundamental  matrix  N  as  depicted  in  Equation  (4),  N  =  (I  -  Q)'1 


-0.9877 

-0.0429 

-0.0413 

-0.491 

-0.0045 

-0.0037 

-0.0131 

-0.1491 

0.1503 

0 

-0.0122 

0.9572 

0 

-0.4904 

-0.0045 

-0.0037 

-0.0131 

-0.1489 

-0.1501 

-0.0424 

-0.0118 

-0.0413 

0.9603 

-0.4727 

-0.0043 

0 

-0.0126 

-0.1435 

-0.1447 

-0.0409 

-0.0118 

-0.0412 

-0.0396 

0.5288 

-0.5288 

-0.0043 

-0.0125 

-0.1431 

-0.1443 

-0.0408 

-0.0137 

-0.048 

-0.0462 

-0.5496 

-0.5496 

0.995 

-0.0146 

0 

-0.1683 

-0.0476 

-0.0118 

-0.0412 

-0.0396 

-0.4712 

-0.4712 

-0.0043 

-0.0125 

-0.1431 

-0.1443 

-0.0408 

-0.0118 

-0.0414 

-0.0398 

-0.474 

-0.474 

-0.0043 

-0.9874 

-0.1439 

-0.1451 

-0.041 

-0.0118 

-0.0412 

-0.0397 

-0.4721 

-0.4721 

-0.0043 

-0.0126 

0.8566 

-0.1445 

-0.0408 

-0.0126 

-0.044 

-0.0423 

-0.5036 

-0.5036 

-0.0046 

-0.0134 

-0.1529 

0.8458 

-0.0436 

-0.0118 

-0.0411 

-0.0396 

-0.471 

-0.471 

-0.0043 

-0.0125 

-0.143 

-0.1442 

0.9592 
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11 

12 

13 

14 

15 

O') 

16 

17 

18 

19 

20 

11 

"1.1515 

0.5302 

0.4882 

6.0698 

0.555 

0.0437 

0.1616 

1.8339 

1.8583 

0.4763' 

12 

0.1515 

1.5302 

0.447 

6.0698 

0.555 

0.0439 

0.1616 

1.8339 

1.8583 

0.5187 

13 

0.1515 

0.5302 

1.4882 

6.0698 

0.555 

0.0402 

0.1616 

1.8339 

1.8583 

0.5187 

14 

0.1516 

0.5304 

0.4882 

7.0722 

0.556 

0.0437 

0.1617 

1.8346 

1.859 

0.5189 

N  =  15(f) 

0.1514 

0.53 

0.4884 

6.0674 

0.555 

0.0437 

0.1615 

1.6663 

1.8576 

0.5185 

16 

0.1516 

0.534 

0.488 

6.0722 

0.556 

1.0437 

0.1617 

1.8346 

1.859 

0.5189 

17 

0.1524 

0.5336 

0.4884 

6.1081 

0.559 

0.044 

0.1626 

1.8455 

1.87 

0.522 

18 

0.1518 

0.5315 

0.4913 

6.0841 

0.557 

0.0438 

0.162 

1.8382 

1.8627 

0.52 

19 

0.162 

0.5669 

0.5219 

6.4894 

0.594 

0.0467 

0.1728 

1.9607 

2.9868 

0.5546 

20 

v0.1515 

0.5302 

0.4882 

6.0698 

0.555 

0.0437 

0.1616 

1.8339 

1.8583 

0.5187 

The  fundamental 

matrix 

N  will  be  interpreted  alongside  the 

associated  variance  b 

(12) 


variance  of  visit  as  in  Equation  (5) 

3.2  Calculation  of  Variances  and  Standard  Deviation 

Variance  on  number  of  visits 
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Nsq  =  (N)2  = 


(14) 
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1.3259  0.2811  0.2383  36.8424  0.0031  0.0019  0.0261  3.3631  3.4533  0.2268 

0.0229  2.3416  0.1998  36.8424  0.0031  0.0019  0.0261  3.3631  3.4533  0.2691 

0.0229  0.2811  2.2147  36.8423  0.0031  0.0016  0.0261  3.3631  3.4533  0.2691 

0.023  0.2814  0.2385  50.0156  0.0031  0.0019  0.0261  3.3658  3.456  0.2693 

_  0.0229  0.2809  0.2381  36.8134  1.1141  0.0019  0.0261  2.7765  3.4506  0.2689 

f<?_  0.023  0.2814  0.2385  36.8712  0.0031  1.0894  0.0261  3.3658  3.456  0.2693 

0.0232  0.2847  0.2143  37.3087  0.0031  0.0019  1.3517  3.4057  3.497  0.2725 

0.0231  0.2825  0.2394  37.0162  0.0031  0.0019  0.0262  8.0554  3.4696  0.2704 

0.0262  0.3213  0.2724  42.1118  0.0035  0.0022  0.0298  3.8442  8.9207  0.3076 

,  0.029  0.2811  0.2383  36.8424  0.0031  0.0019  0.0261  3.3631  3.4533  2.3066 


(16) 


The  variance  of  number  of  visit  N2  =  N  (2Ndg  -  I)  -  Nsq  or  variance  of  average  number  of 
times: 


11 

12 

13 

14 

15 

0) 

16 

17 

18 

19  20 

11 

4).  1744 

0.8114 

0.7265 

42.941 

0.0586 

0.0456 

0.188 

5.2129 

5.789 

0.7435 

12 

0.1744 

0.8114 

0.6836 

42.941 

0.0586 

0.0458 

0.188 

5.2129 

5.789 

0.7878 

13 

0.1744 

0.8114 

0.7265 

42.941 

0.0586 

0.0421 

0.188 

5.2129 

5.789 

0.7878 

14 

0.1745 

0.8116 

07267 

42.943 

0.0587 

0.0457 

0.1881 

5.2136 

5.7899 

0.788 

15h) 

0.1744 

0.8112 

0.7263 

42.939 

0.0586 

0.0456 

0.188 

5.0157 

5.7881 

0.7876 

16 

0.1745 

0.8116 

0.7267 

42.943 

0.0587 

0.0457 

0.1881 

5.2136 

5.7899 

0.788 

17 

0.1754 

0.8147 

0.7296 

42.978 

0.059 

0.0459 

0.1891 

5.2244 

5.8036 

0.7911 

18 

0.1748 

0.8126 

0.7276 

42.955 

0.0588 

0.0457 

0.1884 

5.2172 

5.7945 

0.789 

19 

0.1848 

0.8467 

0.7591 

43.187 

0.0625 

0.0486 

0.1991 

5.3147 

5.934 

0.8224 

20 

.0.1744 

0.8114 

0.7265 

42.941 

0.0586 

0.0456 

0.188 

5.2129 

5.789 

0.7878 

3.3  Interpretation  of  the  fundamental  matrix  N  Eq.(4),  and  its  derivative  N2. 

Essentially,  as  discussed  earlier,  N  estimates  the  number  of  times  or  visit  subjects  (staff) 
starting  from  any  of  the  non-absorbing  state,  transit  to  other  states  within  the  transient  state 
before  being  absorbed,  with  the  associated  variability  estimate  as  given  by  matrix  -  N2  will 
be  interpreted  together. 

Sample  interpretation  from  N  and  N2  as  in  Eq(12)  and  Eq(17)  respectively. 

(i)  20  ->  12  (Near  miss  (Nm)  Minor  Injury  Im) 

N  =  0.5302,  N2  =  0.8114 

The  significance  of  this  entry  is  that  for  every  1,000  near  miss  visit  in  transiting  to  minor 
injury  (Im),  530  times  of  the  near  misses  results  into  minor  injury  (Im)  if  the  safety  function 
remains  constant.  This  can  be  interpreted  also  that  there  is  53%  possibility  of  every  near  miss 
leading  to  a  minor  injury  (Im).  The  event  happens  with  an  associated  variance  of  0.8114  and  a 
standard  deviation  a  =  v/O.8114  =  0.9008  6N2 

(ii)  19  — ►  12  (Unsafe  Act  (UA)  -»•  Minor  Injury  (Im)) 

N  =  0.5669,  N2=  0.8469 

This  entry  signifies  that  for  instance  every  1 ,000  visits  or  number  of  time  of  Unsafe  Act 
transiting  to  minor  Injury,  567  times  will  end  up  in  minor  injury  (Im)  or  57%.  This  is  expected 
because  of  the  volatility  of  the  gases  environment.  This  happens  with  an  associated  variance 
0.8469,  standard  deviation  a  =  V0.8469  =  0.9202  GN2 
Similarly 

(i)  18  —►  17  (N,8,17)  (Unsafe  condition  (UC)  — >  First  Aid  Case(FAC) 

N18>17  =  0.162,  N2(i8,i7)  =  0.1884 
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The  significance  of  this  entry  is  that  in  every  100  times  of  unsafe  condition  exposure,  16  of 
such  exposure  will  lead  to  field  staff  sustaining  injury  that  will  require  first  aid  treatment. 
This  happens  with  an  associated  variance  of  N2  =  0.1884  and  a  standard  deviation  a  = 
VO. 1884  =  0.434  6  N2 

Similarly,  other  essential  inferences  can  be  deduced  from  the  fundamental  matrix  N  and  its 
associated  variance  matrix  N2  for  other  state  parameter  of  interest. 

Computation  and  interpretation  of  the  absorbing  state  matrix 
B  =  (I  -  Q)'1  R  =  NR  as  in  Equation  (8),  N  and  R  matrices  were  stated  earlier. 

12  3  4  5  6789  10 
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^0 

0.005 

0.0027 

0.1161 

0) 

0.0656 

0.7048 

0 

0 

0.026 

0.0798 

2 

0 

0.005 

0.0027 

0.1161 

0.0656 

0.7048 

0 

0 

0.026 

0.0798 

3 

0 

0.005 

0.0027 

0.1161 

0.0656 

0.7048 

0 

0 

0.026 

0.0798 

4 

0 

0.005 

0.0023 

0.1162 

0.0657 

0.705 

0 

0 

0.026 

0.0798 

NR  =  5  (i) 

0 

0.005 

0.0027 

0.1161 

0.0656 

0.7045 

0 

0 

0.0264 

0.0797 

6 

0 

0.0046 

0.0027 

0.1162 

0.0657 

0.705 

0 

0 

0.026 

0.0798 

7 

0 

0.0051 

0.0027 

0.1169 

0.0661 

0.7092 

0 

0 

0.0261 

0.074 

8 

0 

0.005 

0.0027 

0.1164 

0.0658 

0.7064 

0 

0 

0.0237 

0.08 

9 

0 

0.0054 

0.0029 

0.1242 

0.0702 

0.6843 

0 

0 

0.0278 

0.0853 

10 

1° 

0.005 

0.0027 

0.1161 

0.0656 

0.7048 

0 

0 

0.026 

0.0798 

f 


B  = 


Fatality 

Permanent  Disability 
Injury  (Severe) 


—>  0  Third  Party  Fatality  -»  0.005 N 

0.0027  Fire  Incidence  (Severe)  0.1161 

-»  0.0656  Medical  TreatmentCase  (Severe)  -»  0.7048 


[19] 


RoadTrafic Accident (Severdd  0  Oil  Spill  (Call Outs)  0 

VR estrictedWork  Cases  (Severe)  -»  0.026  HumanError  -»  0.0798y 

The  numbering  of  matrices  N,  N2  and  B  is  in  consonant  with  the  states  as  depicted  in  Table  1 
which  facilitates  cross  references.  The  B  matrix  gives  the  long-run  transition  of  subject  (staff) 
within  the  system  which  shows  the  general  trend.  It  is  evident  that  the  row  entries  in  each 
column  are  the  same  showing  that  it  is  actually  represents  a  stabilized  matrix  (trend)  and  thus, 
can  be  interpreted  column  wise.  Column  1  represents  state  1  which  is  fatality,  while  column  2 
represents  state  2  which  is  third  party  fatality,  just  to  mention  a  few.  For  example,  column 
wise,  column  6  which  represents  medical  treatment  case  Severe  (MTCs)  reveals  that  704 
subjects  in  every  1,000  are  going  to  be  categorized  to  have  a  severe  medical  condition  which 
requires  treatment  i.e.  (medical  treatment  case  severe)  if  the  trend  of  injury,  illness  and 
gaseous  exposure  remain  the  same.  While  fire  incidence  severe  in  column  4  implies  that  in 
every  100  subject,  about  11.6  i.e.  12  will  lead  to  severe  fire  incidence  (  FIS). Plausible  too! 

The  accompanying  graphical  representation  in  matrix  Equation  (18)  maps  states  to  the  long- 
run  probabilities  of  occurrence  as  discernable  from  the  B-matrix  depicted  supra. 


The  total  movements  of  subject  within  the  non-absorbing  state  before  being  absorbed  is 
computed  as  follow.  This  movement  is  represented  by  the  r-matrix  given  by 
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v0.152 

0.53 

0.49 

6.07 
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1.83 

1.858 

1.519  ) 

X 

v12.712. 

And  its  associated  variance  is  T2  as  in  Eq  (7)  and  matrix  Equation  (22) 
r2=  (2N  - 1)  t  -  rsq 
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The  Eq.(20)  for  r  suggests  that  subjects  on  the  average  change  position  or  habituate  about 
13times  among  the  non-absorbing  states  before  being  finally  absorbed  (trapped)  into  any  of 
the  ten  absorbing  state  numbering  1-10  in  Table  1.  While  the  accompanying  matrix  (xxii) 
estimates  the  associated  variance  expected  number  of  movement  or  habituation  T2  =  152  and 
standard  deviation  Eq  (23)  of  about  12.  This  implies  that  the  estimation  of  N  and  B  computed 
could  hover  about  this  mean  value  of  12. 

Transient  probability  of  visiting  transient  states  refer  to  Equation  (9)  H  =  (N  - 1)  Ndg'1 


Transient  probabilities  of  visiting  transient  states  H  =  {N  -l)Ndg  1 
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(25) 


The  H-  matrix  as  in  matrix  Equation  (25)  estimates  the  probability  (chances)  of  a  subject 
transiting  amongst  the  non-absorbing  state.  For  instance  H2o,  i4  =  0.8583  the  transient  to 
transient  estimate  or  likelihood  probability  of  visit  from  near  misses  to  medical  treatment 
case  minor  (MTCm)  has  85  chances  of  occurrence  in  every  100. 

Furthermore,  H 1 9_  19  =  0.665,  the  transient  to  transient  probability  visit  to  unsafe  acts  give  67% 
chance  or  likelihood  probability  which  agrees  with  the  high  state  value  total  in  Table  1 

4.  Conclusion 

Safety  performance  function  is  vital  to  the  growth  and  survival  of  all  engineering  firm  in  term 
of  increase  productivity,  cutting  edge  competitive  advantage  and  over  all  safety  of  the 
employee  and  the  environment.  This  can  be  achieved  by  applying  better  approaches  and 
model  in  predicting  accident  outcomes  which  the  Markov  chain  model  comes  handy. 

Results  from  the  study  suggest  certain  trend  and  pattern  as  depicted  in  the  long-run  transition 
distribution  safety  matrix  B  =  NR  as  subjects  transit  or  habituate  from  a  non-absorbing  state 
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before  entering  any  of  the  absorbing  states.  Fatality,  permanent  disability,  severe  injury  was 
significantly  low.  Although,  there  are  other  absorbing  state  parameters  that  are  of  concern  as 
evident  in  the  study,  such  as  human  error  (HE)  8%  with  an  attendant  severe  fire  incident 
occurrence  of  11.6%  because  of  volatility  (gas)  of  the  working  environment,  which  signals  a 
worrisome  development.  Remarkably,  70.5%  of  the  field  workers  had  severe  medical 
treatment  case  (MTCs)  over  a  16-year  period  traceable  to  occupational  injury  and  illness  due 
to  poisonous  inhalation  of  gases. 

Consequently,  the  significance  of  the  above  import  of  severe  medical  treatment  case  to  the 
organization  under  study  is  weighty  in  terms  of  cost  of  medical  bill,  reduced  productivity,  lost 
time  days  and  increase  workload  for  staff  replacing  the  affected  co-worker.  From  the 
foregoing  managers  can  gain  deeper  insight  of  incident  level  from  the  application  of  the 
Markov  chain  model  as  a  veritable  tool  for  decision  making.  For  instance,  the  transient  to 
transient  movement  (H)  and  other  transitions  or  habituation  to  absorbing  states  made  by 
workers  can  be  used  to  signal  impending  danger  or  can  be  shown  on  bill  board  as  caution  to 
workers  on  how  they  transit  or  habituate  in  the  various  states  before  entrapped  in  an 
absorbing  state,  with  the  right  remediation  in  place. 

In  summary,  the  Markov  chain  model  offers  a  veritable  "framework  for  the  development  of 
safety  system  for  similar  organization. 
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